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A number of the proposals for supporting a Galileo S-band (2.3-GHz) mission 
involve arraying several antennas to maximize the signal-to-noise ratio (and bit rate) 
obtainable from a given set of antennas. Arraying is no longer a new idea, having 
been used successfully during the Voyager encounters with Uranus and Neptune. 
However, arraying for Galileo’s tour of Jupiter is complicated by Jupiter s strong 
radio emission, which produces correlated noise effects. This article discusses the 
general problem of correlated noise due to a planet, or other radio source, and applies 
the results to the specific case of an array of antennas at the DSN’s Tidbinbilla, 
Australia, complex (DSS 42, DSS 43, DSS 45, and the yet-to-be built DSS 34). 
The effects of correlated noise are highly dependent on the specific geometry of 
the array and on the spacecraft-planet configuration; in some cases, correlated 
noise effects produce an enhancement, rather than a degradation, of the signal- 
to-noise ratio. For the case considered here— an array of the DSN’s Australian 
antennas observing Galileo and Jupiter— there are three regimes of interest. If the 
spacecraft-planet separation is <75 arcsec, the average effect of correlated noise is 
a loss of signal to noise (~0.2 dB as the spacecraft-planet separation approaches 
zero). For spacecraft-planet separations >75 arcsec, but <400 arcsec, the effects 
of correlated noise cause signal-to-noise variations as large as several tenths of a 
decibel over time scales of hours or changes in spacecraft-planet separation of tens 
of arcseconds; however, on average its effects are small (<0.01 dB). When the 
spacecraft is more than 400 arcsec from Jupiter (as is the case for about half of 
GaIileo } s tour), correlated noise is a <0.05-dB effect. 


I. Introduction 

It is now becoming common practice to array a number 
of antennas for telemetry reception in order to increase the 
effective aperture [1,2]. The signal-to-noise ratio (SNR) 
improvement obtained by arraying is straightforward to 


calculate when a pointlike spacecraft is the only source in 
the array’s field of view. However, the calculation becomes 
complicated when an additional hot body (such as a radio- 
loud planet) is in the beam since, unlike system noise, 
the incoming radiation from this hot body is correlated 
at different antennas. This article presents a calculation 
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of the effect of correlated noise on the SNR achieved by a 
given array configuration. Other aspects of this topic have 
been discussed in [3] and in references therein. 

This article is organized as follows: Sections II and III 
present a general formulation of the problem, Section IV 
considers the special case of an array of N identical anten- 
nas, and Section V looks at the specific case of arraying the 
DSN complex in Australia for Galileo S-band (2.3-GHz) 
telemetry reception. To make the article more readable, 
many of the detailed calculations have been relegated to 
the appendices. Much of the discussion parallels Hjellm- 
ing’s treatment of the Very Large Array (VLA) [4]. 


II. Outline of the Problem 

Consider an array of N (not necessarily identical) an- 
tennas. The on-axis gain 1 of the ith antenna is denoted by 
Gi and its system temperature (looking at an empty field) 
is denoted by 7}. The direction toward a point on the sky 
can be described by a unit vector s, with Sq pointing to- 
ward a source at the array's field center. 2 The shape of the 
tth antenna beam can be described by the field pattern, 3 
a dimensionless function /*( s - s 0 ) such that f?( s - s 0 )Gi 
describes the effective gain of the ith antenna in a direc- 
tion s — s 0 from the beam center. In general, G,-, 7}, and 
/• may all be elevation dependent. This dependence is not 
considered here; it can easily be included in quantitative 
calculations. The geometry of the array of N antennas 
can be described by the set of baseline vectors between 
each antenna pair Bj* = — B*,*, where B,* denotes the 
vector from the ith antenna to the *th antenna. The basic 
geometry is shown in Fig. 1. 

It is assumed that the spacecraft is transmitting a 
narrow-band, polarized signal with a total Earth-received 
power per unit area of P s/c and that the planet has a to- 
tal, unpolarized flux density (power per unit area per unit 
bandwidth) S>. The angular distribution of the planet’s 
radio emission is described by the brightness distribution 

1 Throughout this article, an astronomer’s definition of antenna gain 
is used (commonly measured in units of K per Jy) of G = zA/2k By 
where c is the dimensionless antenna efficiency, A is the physical 
antenna area, and k B is Boltzmann’s constant. 

2 Here, i 0 is defined as the nominal direction in which the antennas 
in the array are pointing, which is the direction assumed when 
model delays are calculated (see Appendix A). As discussed in 
Appendix C, it is not necessarily the phase center of the array. 

In general, the field pattern is considered to be a complex quantity 
(e.g., [5, p. 27]), but for the purposes of this article, it is sufficient 

to assume that it is real. 


(flux density per unit solid angle) 7/>(s), which is assumed 
to be constant in time and only slowly varying with fre- 
quency (i.e., constant across the observing bandwidth). 

Fairly simple assumptions are made about the method 
of arraying: The incoming RF signal at the ith antenna is 
mixed to baseband, after which a phase shift is inserted. 
This phase shift has a component <£[" (the model phase) 
that can be used to compensate, in real time, for phase 
differences between antennas (see Appendices A and C for 
details). The baseband signal is delayed by a model delay 

, and finally, the delayed signals from the N antennas 
are summed with appropriate weighting (see Appendix D). 
It is assumed that the signal is filtered through a rectan- 
gular baseband filter of width Av; for an upper sideband 
system (which is assumed throughout), the sky frequency 
at the center of the passband is v\ ky = i/ io + Av/2, where 
U]° is the mixer frequency. It is also assumed that the 
model delays are chosen so as to compensate for the delay 
that a signal from the direction sq (the field center) suffers 
between an arbitrary reference point and the ith antenna. 



where r, is the vector from the reference point to antenna 
i. A simple block diagram of the signal path is shown in 
Fig. 2. 

The total arrayed power can be written as the sum of 
the arrayed power from each of the three sources P fc = 
"F > where p£ c } P £ , P£ , represent, res pec- 

tively, the arrayed power from the spacecraft, from the 
planet, and from system noise. 4 The effective SNR of the 
array is proportional to the ratio of the arrayed space- 
craft power to the total arrayed power; in the low signal 
limit, this reduces to the ratio of P^ c to the power re- 
ceived from all sources other than the spacecraft (in this 
case, the planet and system noise). A figure of merit f3 can 
be defined for an array; f3 is analogous to G/T for a sin- 
gle dish and quantifies how, for a given spacecraft power 
(' p ’/c ) an observing bandwidth (At/), the SNR varies with 
the following array configuration: 

Ai/ P^ c At/ Pi! c 
W,/c Ps * 2V,, C Pi + Pi (2) 

4 System noise includes receiver noise, sky background, and ground 
pickup— all sources of unwanted noise other than individual, iden- 
tifiable sources (such as a planet) in the field of view. 
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The term P^'/Ps is the SNR; the term Au/2 P’J is a 
normalization factor that removes from /? any dependence 
on quantities other than array and source geometry. 


III. General Expression for p 

As discussed in Appendix A, the summed baseband 
voltage can be written in the form 


spacecraft (P t >,c , p‘1% the planet {Pf , p£ k ), and random 
noise (P? , p? k ), so Eq. (4) can be written 


Pz=jr w W e + p ? +p '” ) 

« = 1 

+<£+<>!!> « 

t=l k?i 


= ^2 Wie -1 *" Vi(v,t + if) (3) 

« = 1 

where tf) is the baseband voltage given by Eq. (A-4), u 
is the baseband frequency, and Wi is the weighting factor. 

The arrayed power is then 

Ps = (^J 


Thus, P E , the total power in the arrayed signal, can 
be written as the sum of contributions from the three 
sources (the spacecraft, the planet, and system noise), 
p s = p' /c + P£ + with 

(8) 

i - 1 »=1 

p£ = J2 W?P? +Y.Y. WiW k rf k (9) 

i=l i= 1 





AT N N 

= Y^wiPi+Y1^2 WiWkPik 

»=1 *=1 M* 

where the angle brackets denote time averaging and T is 
the interval over which the time average is taken; P, is the 
single dish power from the ith antenna, 

Pi-lj du |u,-(i',<)| 2 ^ (5) 

and pik is the unnormalized correlation between the ith 
and Arth antennas, 


Pg + t t WiW kP ? k 

i=l 1 = 1 k*i 

( 10 ) 

«=i 

Because the phases of the individual p£ k terms in Eq. (10) 
are random, their sum is always small (on average, zero) 
and was dropped. However, the sums over p£ k and p ik are 
not necessarily small; in fact, a key aspect of arraying is 
choosing the inserted model phases to ensure that the sum 
over p\ k c is maximized (see Appendix C). Expressions for 
p.*/ c pP t and P-* are derived in Appendix A, Section II, 
and given’ by Eqs. (A-13), (A- 14), and (A-15), respectively. 
Expressions for p’^ and p£ k are derived in Appendix A, 
Section C, and given by Eqs. (A-25) and (A-28). Substi- 
tuting Eqs. (8)— ( 10) into Eq. (2) gives 


jrw?p; /c + 

A 1/ »•=! » = 1 

2V~T X ~n n n 

J2 w ? ( p i N + P > p ) + £ £ WiW kt f k 

i = l « = 1 

As can be seen in Eqs. (A-ll) and (A-18), both Pi and 

p ik can be written as the sum of contributions from the '■ 


Pik = (J o V + rf* )«:(!/,* + rf)) a 


( 6 ) 


0 = 
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The phases of p*J k c and pf k depend on the quantity 6<f> { — 
S<t>k — — &T “ <f>k 4- <t>™, where are the actual 

antenna-based phases, which include both hardware and 
atmospheric effects, and are the model phases. 

If no attempt is made to compensate for the individual 
antenna phases (i.e., if ^ = </>% = 0), the phases of p*[ c 
and pf k are uncorrelated from baseline to baseline, and in 
both Eqs. (8) and (9) the sums over the cross-correlation 
will be small (zero on average). However, as discussed in 
Appendix C, the total phase of p s J k c will be zero for all i, k , 
if the values of (fr™ ,4>™ are chosen to satisfy Eq. (C-l): 


can become significantly nonzero. This is the source of the 
correlated noise contribution from the planet. 

Equations (A-13), (A-14), (A-15), (A-25), and (A-28) 
can be used to expand Eqs. (8), (9), and (10), and if the 
array is phased on the spacecraft [Eq. (C-l) is satisfied], 
the arrayed fiower from the spacecraft, the planet, and 
system noise can be written, respectively, as 


P^ = 2k B V, /c 


N N N 


1=1 t=l (fcjfi 


( 12 ) 


Hi - Hk = 4>i - 4>? - 4>k + <t>? 

= -2„S/'*“.[s, /e -s 0 ] 
C 


P£<t> = Ai/kpSp 


N 


E 


When are chosen to satisfy the above expression, 

all the terms in the double sum in Eq. (8) add in phase, 
and for large N , the sum over p\ k c contributes signifi- 
cantly more t o P^J C than does the sum over P*J C (since 
pik ~ y/PiPjc). Phases chosen to satisfy Eq. (C-l) max- 
imize P ^ c , and except in pathological 5 cases this choice 
maximizes /?. The process of determining and inserting 
these values of 0™ is referred to as phasing the ar- 
ray. In the following, the subscript <f> is used to denote 
expressions that assume that the array is phased on the 
spacecraft. 

Unfortunately, the values of <f > which maximize 
the sum over p*£ c y may also de-randomize the phases of 
as a consequence, sums over pf k in Eqs. (4) and (11) 


N N 

+ £ £ Wi W k f ip h P Tf k 

i=l k?i 


(13) 


N 

P$* = ( 14 ) 

1=1 


In the above expressions, fi P , h P represent the average 
value of the beam pattern in the direction of the planet, 
and Tf k [defined by Eq. (A-27)] is a dimensionless, complex 
quantity, with magnitude less than or equal to unity, which 
depends only on the array geometry and the structure of 
the planet (or other background source). 

By inserting Eqs. (12), (13), and (14) into Eq. (2), the 
expression for /? reduces to 



N N N 

£ W?Gi : + £ £ W < Wk 
/=! i=l k& 


Sp 

> =1 

N N 

Gi + £ £ WiW k f ip f kp y/G&rr k 
*= 1 *** 

i=l 


5 One can imagine cases where the increase in P^J C obtained by 
phasing the array is more than counteracted by the resulting in- 
crease in P £ , but even attempting to array in such cases would be 
difficult. 


IV. An Array of Identical Antennas 

The major effects of correlated noise are easily seen by 
examining an array of identical antennas. Consider an ar- 
ray of N antennas, each with a gain G, a system tempera- 
ture T, and an average field pattern in the direction of the 
planet fp . Since the antennas are identical, Wi = 1 for all 
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( 22 ) 


*. When this array is properly phased on the spacecraft, 
Eqs. (12), (13), and (14) become 




NG 

fpGSp + T 


p£ = 2 k B N 2 GV,, c 
P£<t> = Avks/pG Sp 


N N 
* = 1 k^i 


Pvt = Auk B NT 


and Eq. (15) becomes 


NG 



N N 


f 2 P GSp 

i + n-'Z 

+ T 




(16) 

(17) 

(18) 


(19) 


Since \P[ k \ < 1 for all baselines, it is always true that 

i+ r ‘EE f ^ w (20) 

* = 1 kjU 


This is equivalent to observing the spacecraft with a single 
dish of gain NG and an effective system temperature of 
P p GSp+T. 

In both cases, the effective gain is NG ; it is indepen- 
dent of baseline length and depends only on the gains of 
the individual antennas. However, the effective system 
temperature is not the same in the two extreme configu- 
rations. The planet contributes N fpGSp to the system 
temperature in the compact array limit but only fpGSp 
in the extended array limit. In the extended array limit, 
the noise contribution from the planet is simply the sum 
of its (uncorrelated) contributions to the individual system 
temperature. In the compact array limit, its contribution 
is N times larger due to the correlated noise effects. 

Though, as discussed below, the effects of correlated 
noise are a complicated function of array geometry. In gen- 
eral, the more extended the array configuration , the smaller 
the correlated noise contribution of a planet or other hot 
body. As a rule of thumb, correlated-noise contributions 
are significant on baselines where ‘'tk y BikRp/ c £ 1 (see 
Appendix A, Section III). 


From Eq. (A-27) [see also Eqs. (B-l), (B-3), and (B-5)], 
it can be seen that, in the limit of an extremely compact 
array [i.e., Bik << c /i l/ sky^P^’ f° r *>^> being the 
characteristic angular radius of the planet], T P k — * 1 on all 
baselines and 


It should be noted that in the intermediate cases where 
Bik ~ C/(v c sky R P ), the sum over T? k may, for some ge- 
ometries, be negative. In such cases, the performance of 
the array is actually enhanced over that of the extended 
array limit. 


NGSp + T * 

where the substitution fp = 1, appropriate in this limit, 
has been made. This is equivalent to observing the space- 
craft with a single dish of gain NG and an effective system 
temperature of NGSp + T. Not surprisingly, observing the 
spacecraft and planet with a compact array is analogous to 
observing the two objects with an antenna with N times 
the area of a single array element. 

In the limit where B ik » c/(v c sky R P ), for all i,k (i.e., 
an extremely extended array), \Ffk\ — ► 0 on all baselines 
and 

6 If the baseline is small compared to c/j/^Kp), the diameters of 
the individual antennas would be smaller still. 


V. DSN Complexes at S-Band 

One of the proposals for support of the Galileo S-band 
mission is the arraying of antennas at each DSN com- 
plex. This section describes the performance of an array 
at the Canberra complex that includes DSS 42, 43, 45, 
and the soon-to-be-built DSS 34. Table 1 lists each an- 
tenna, its S-band gain and system temperature at zenith, 
its S-band beamwidth at full-width half-power (FWHP), 
and its station coordinates (east, north, and vertical) rel- 
ative to DSS 43. 

With the model of Jupiter given in Appendix B, Eq. (15) 
can be used to calculate /fy for the array. The improvement 
that the array would provide relative to the stand-alone 
use of DSS 43, if correlated noise effects could be neglected, 
is given by the ratio 
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N N N 

Y w ?Gi + YH WiWk \/^r 

i = 1 1 = 1 /b^i 


N 

I 

i = l 


AT 


5p E ^. 2 /.p g - + E W ? T < 


i=i 


^43 + /Isp^a^p 
G 43 


(23) 


aration between the spacecraft and the center of Jupiter 
are held fixed. Each point on the plot then represents 
the value of (3$/ 0$^ for a randomly chosen value of the 
hour angle, the orientation on the sky of Jupiter's radi- 
ation belts, and the orientation of the spacecraft-Jupiter 
separation. Thus, the density of points on a portion of a 
plot provides an estimate of the likelihood, as a function of 
the hour angle, of obtaining a particular value of 0$ / 0^ 
for the given parameters. 


When correlated noise is neglected, this ratio is indepen- 
dent of hour angle, array geometry, or source structure, 
but it depends on Sp — the flux of Jupiter — and, there- 
fore, on the Earth-Jupiter distance. At opposition, this 
distance is 4.2 AU, the S-band flux of Jupiter is approxi- 
mately 5.8 Jy and, with the parameters given in Table 1, 
0^ / /?43 = 1.41. At conjunction (an Earth-Jupiter dis- 
tance of 6.2 AU), Sp « 2.6 Jy and 0^/043 — 1.38. The 
improvement is not as great at large Earth-Jupiter dis- 
tances because Jupiter's contribution to the system tem- 
perature is less; for comparison, if Sp = 0, /?<£«, / 043 = 
1.34. It is useful to note that despite the effects of corre- 
lated noise, arraying is most useful when extended back- 
ground sources are present, particularly if the baselines 
are long. In the above calculations, it is assumed that 
fp = 1 for sill antennas and, as in all the calculations in 
this section, W> = y/Ui/Ti [see Eq. (D-4)]. Throughout 
this section, it is also assumed that the gains G, and sys- 
tem temperatures 7} do not vary with elevation, which 
leads to a slight overestimate of correlated noise effects at 
low elevations. 

To assess the effects of correlated noise, the SNR pro- 
vided by the actual array is compared with that provided 
by the same antennas if arrayed in a configuration with in- 
finitely long baselines, examining the ratio 04 >/ 0 ( p oo . This 
ratio depends not only on Earth-Jupiter distance, but, in 
a complicated way, on the array and source geometries 
and on hour angle. Figure 3 plots this ratio for a num- 
ber of different geometries. In each plot the declination of 
Jupiter, 7 the Earth-Jupiter distance, and the angular sep- 


7 A declination of —21 deg, corresponding to that of the Galileo en- 
counter in December 1995, has been used for Fig. 3. The results for 
other Jupiter declinations are different in detail but qualitatively 
very similar. 


Figures 4(a) and (b) summarize the results of the 
calculations shown in Fig. 3 plotting, as a function of 
spacecraft-Jupiter separation, the average (over hour an- 
gle and orientation) value of 0$ / 0 $^ , as well as its min- 
imum and maximum values. Thus, each plot in Fig. 3 is 
reduced to three points in Fig. 4 — one on an average curve, 
one on a maximum curve, and one on a minimum curve. 
The following conclusions can be drawn from Fig. 4: 

(1) When Galileo is within ~75 arcsec of Jupiter, the 
average effect of correlated noise is a loss of SNR, 
which may be as large as 5 percent at zero separa- 
tion. 

(2) When Galileo is more than 400 arcsec from Jupiter, 
the effects of correlated noise are negligible. 

(3) For separations in the range of 75 to 400 arcsec, the 
effects of correlated noise are small on average. How- 
ever, the loss may be significant for certain geome- 
tries (as may be the enhancement). In this range of 
planet-spacecraft separations, careful calculations of 
correlated noise are necessary in a case where a 5- 
percent SNR loss would be critical. 

The above conclusions refer specifically to an array con- 
sisting of the antennas listed in Table 1. Correlated noise 
effects are likely to be similar for arrays of similar antennas 
with similar baseline geometries, but quantitative calcula- 
tions for other arrays have not been carried out. Because 
correlated noise causes significant SNR variations, and be- 
cause these variations are highly geometry dependent, de- 
tailed calculations of these effects should be done for any 
situations in which a few tenths of a decibel of SNR are 
crucial. For most purposes involving Galileo and Jupiter, 
Eqs. (15), (B-6), and (B-9) should be suitable for such 
calculations. 
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Table 1. S-band parameters of antennas In proposed Canberra array. 


Antenna 

Gain, 

K/Jy 

System 

temperature, K 

S-band beamwidth 
(FWHM), deg 

Station coordinates 1 

Bast, m 

North, m 

Vertical, m 

DSS 43 

0.95 

18.5 

0.11 

0.0 

0.0 

0.0 

DSS 42 

0.21 

22.0 

0.28 

0.0003 

194.1921 

-13.6414 

DSS 45 

0.16 

38.0 

0.23 

-325.3907 

440.1822 

-13.1378 

DSS 34 b 

0.16 

30.0 

0.23 

68.8 

440.2 

0.0 


1 Relative to DSS 43. 
b Values for DSS 34 are approximate. 
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ARRAY 

REFERENCE 

POINT 


Fig. 1. The basic geometry of an array. Antennas / and k are lo- 
cated at, respectively, r/ and r* with respect to the array reference 
point and are pointed toward the direction $o. It is clear from the 
diagram that the delay suffered by a wavefront from a direction £ 
between antenna / and antenna k is proportional to Bfo • $. 
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Fig. 2. A simplified block diagram of an array signal path. 













HOUR ANGLE, hr 


Fig. 3. The ratio of ( 3 $/ ( 3 ^ for an array of antennas at the DSN Canberra complex observing Jupiter and a spacecraft at S-band, 
for a variety of array and source geometries. The declination of Jupiter = -21.0 deg and the array consists of DSS 43, DSS 42, 
DSS 45, and DSS 34. 
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(a) 



JUPITER-SPACECRAFT SEPARATION, arcsec 

Fig. 4. The average, minimum, and maximum values of 
obtained as source orientation is varied, plotted as a function of 
the separation between the spacecraft and the center of Jupiter, 
with Jupiter declination = -21.0 deg: (a) Earth-Jupiter distance = 
4.2 AU and (b) Earth-Jupiter distance = 6.2 AU. 





Appendix A 
Basic Expressions 


I. Incoming Signals 

The voltage (as a function of time t and sky frequency u, ky ) induced in a single polarization channel at r,, the focal 
point of the *th antenna, by a distant 8 source (e.g., a planet or spacecraft) can be written 


V i (u tky ,t)=e i * i y/k^G i JJ ds /,(s - s 0 )Z(i/,i yi f,s)e*( 2TI W< - (»**/c)] + 0(u lky , s)) 


(A-l) 


In this expression, ks represents Boltzmann’s constant and c is the speed of light; 6', is the on-axis antenna gain, /, 
is the magnitude of the antenna’s field pattern, s is a unit vector in the direction of the source, l(v, ky ,t, s) is the 
amplitude of the electric field due to the source, 6(i/ aky ,s) is a phase term (which for most astronomical sources can be 
assumed to be uncorrelated over v aky and s), and <6, is an antenna-based phase shift (including atmospheric effects). 

The voltage induced by system noise can be written 

V t N (v, ky ,<) = VhPFie' ( 27r "' fc y < + e ? Kiy-0) ( A -2) 


where 7} is the noise temperature and 6^{u shy ,t) is a random phase that is uncorrelated over time and frequency 
intervals satisfying Ai/ s k y At > 1. The total RF voltage can be written as the sum of terms due to the spacecraft, 
planet, and system noise: 


Vi(v.ky, t) = V‘ /C (v aky ,t) + K' P (^tv,<) + VWt,,*) 


(A-3) 


When the RF signal is mixed to baseband, it is phase shifted by 2irv to T?'; this stops the fringes, allowing time averaging 
over a longer interval. An additional phase shift - <j> J" can be inserted to (attempt to) compensate for the antenna-based 
phase shifts fc. Without this additional phase shift, the baseband voltage can be written 

= e i2 ^r-t] [v;i\ v , ky ,t) + v p 0 ,, ky ,t) + v/> JJtvI t)] 

= JJ ds /,( s - s 0 )e <2 "'*‘»f‘-( i '•/=)] Jj >/c (^ ty)<) s) e **' / ‘ 

+ Ip^.ky , t, i)e ,#P ( , '-*>'' i >J + e i2x Vlo tr,--0 (v.kv.t)) 


{V'k y ,*) 


= JJdi fi(s - s 0 )e i2T (^-< ir -/ c )] +, ''«[ r . m -( ir -/ c )J) \j a , c {v aky ,t, s)e i0,/ ' 




+ lp{^ky,t,sy eP(v ^ i) + 


(A-4) 


8 Here, “distant” means a source sufficiently far away that, over the extent of the array, a plane wave approximation is valid. 
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where v = u >ky - v lo is the baseband frequency. Like the RF voltage, the baseband voltage can be written as the sum 
of terms due to the spacecraft, planet, and system noise: 

Vi(v,t) = vf c {v, t) + vf (u,t) + vf ( v,t ) (A-5) 

If the additional phase shift is inserted, the baseband signal has the form e - ’*- v, •(*/,<). 

II. Single-Dish Power 

For the tth antenna, the single-dish power can be written 



+ v‘ /e (v,t)vf*(v,t) + vf (v,t)vf*{v,t) + vf{v,t)vf*{v,t) 


+ vf c (v,t)vf*{i/,t) + vf (i/,t)vf’(ts,t) + vf (v,i)vi *(^,<) 


(A-6) 


J / r 


where the angle brackets denote a time average and T is the interval over which the time average is taken. Looking at 
each term separately and considering the planet term first, 


*(!/+ Z/ |o( S) 


vf{u,t)vf'{v,t) = e^VklGi j J ds />'(* — r '^ c ^ +1/, ° (r,m-( * ri/c ^^2'p( l/ + t'lo, s)e ifff 

x e-^'Vk^Gi J I ds' fi(s' - s 0 )e" i2 " (l ' [t - (i ' r>/c)1+ ‘ ,, ‘ >tT - m - (i ' r - /,:)]) Ip(^ + «W, s')e-*'>+‘' 1 ‘>'*') 

= k B G { If di Ms - so )l P {u + u lo ,t,s) I f ds' 

^[^(jz+Olo, - S 0 )Ip(l^ + via, t, s')} 


,+^io[8-*8 / ]r l )/c) 


x e 1 


= kgGi JJ ds /^(s — &o)T^{y ~h Wo , t, s) 

= ksGi j j ds f?(s — io) fr (, + ^ s) 


(A-7) 
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where the integral over s' is nonzero only when s' = s because the radiation from the planet is spatially incoherent, 
and where the substitution l\{v -f v Xo ,t, s) = Ip(v + v\ c , t, s) was made; Ip is the brightness distribution (flux per unit 
solid angle) of the planet. Similarly 

v i h (y^) v i /e *iy, 0 = k B Gi JJ ds /,?( 8 - s 0 )I, /c (v+ u, 0 ,t,s) (A-8) 

with I s f c being the brightness distribution of the spacecraft (usually assumed to be pointlike), and 


v?(t/ t t)vV'*(v,t) = k B Ti 


(A-9) 


Compared to these terms, the cross terms are all small, as can be seen by considering 

t '* /e (M)«f ’(M) - k B Gi JJ* ds fi(s-so)l, /c (u + u, 0 ,t,s) JJ ds' 




« 0 (A-10) 

where the integral over s' is small because 8*t c and 9 P are uncorrelated [unlike in Eqs. (A-7) and (A-8), this integral 
is, on average, zero, even when s' = s]. 

One can therefore write 


■-(i: 


dv |v<(^,<)| 2 


f Av 


*'K /e (*',0l a + h>fM)| 2 + | 


«?W)I 2 ^ 


= p’ ,c + p, p + pi 


N 


(A-ll) 


It is assumed that the spacecraft is a point source located at s,/ c , transmitting a narrowband, polarized signal at a 

sky frequency v‘J k e y with a total power V,/ c . The brightness distribution seen by a receiver matched to the signal’s 
polarization 9 is 


h/ci^aky, t > s) — 2V,/ C 6 (^aky ~ U aky) ^(® — ®a/c) 


(A-12) 


(with 6 representing a Dirac delta function), so 


1 The factor of 2 in Eq. (A-12) arises from the assumption of matched polarizations. 
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P t ,fe = 2k B f?(s, /c - so )GiV,, c « IksGiV./c 


(A-13) 


where it is assumed that /?( s a / c — So) « 1 (i e., 


the spacecraft is close to center of the array’s field of view). 


It is assumed that the planet is an extended, unpolarized, broadband source 
of Sp, which is constant in time and varies only slowly with frequency (i.e., 
observing bandwidth Ai/), so /p(^j* y ,<,s) = Ip( s) and 


of total flux at the frequencies of interest 
it can be considered constant across the 


Pf = k B GiAi/ 



ds f? p { s - so)Ip(s) = k B fi p GiAuSp 


(A- 14) 


where /? is a weighted average of ff p over the planet. 
Finally, from Eq. (A-9) 


P? = k B TiAu 


(A-15) 


III. Cross-Correlation 

The cross-correlation of the voltages from the ith and kth antennas is given by 

p ilc = + r^vK^t + T T)) r (A -16 ) 

By using the general form of the baseband signal [Eq. (A-4)] and because the signal is spatially incoherent, this becomes 


Pik 


= / dvjW-*?) [< /c (t /,< + rf) + vf{v,t + rD + vi r (v,t + 7f )] 
x [vl’%, t + rD + t + rD + v* t N (v, t + tJ* )] ^ 


/•Ai' 


dv e 


*0£ 


-♦D + rr"K ,/C (^.< + T k) + »/>■* + rrK ,e (u,t + r?) 


+ + rr)vl 3/c (v,t + rD + v’ ,c (v,t + Tf>; P (".* + r D 

+ »/>,< + r/ n K P (t',< + r fc m ) + v?(v,t + T?)vl p {v,t + rD 

+ v 3/c (p,t + r, m )vl N (p,t + rD + v? (i /,< + T?)vl N {y,t + r fc m ) 


+ + r in] y 


(A- 17) 
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This can be written as 


where 


Pik = Pik +Pik+ Pik 


(A- 18) 


sfc 

Pik ~ 


P?k = 




= ^ jf dvt$*k-*~) v ;/ e ( Vi t + rf K’ /e (i/,< + rf 

ft = ( Jo dueKK ~ K)v i ("-* + ^)V(",' + T k'Jj 

*?(•'>* + r />i JV (". < + td + «/ w + rn®; ,/e (^, < + td 


(A-19) 


(A-20) 


^eW-C) 


+ u/V,i + »TK ,/, (v,( + rf) + < /e (t/,< + TfK'V,* + rf) 

+ «,*(!/, t + r/ n ) V ; /> (i/, t + rf) + v' /c (i/,< + rfK'V,* + r t m ) 

+ t;f( i /,< + rrK /V (^< + r i m )]\ 


(A-21) 


These expressions are similar in structure to the expression for single-dish power [Eq. (A-6)]. Looking at pf k , one finds 
”/>.* + rrH P (u,t+rn = e^’-^keVG^e^+^r-^) ^JJ </i | e -(2'l*'+‘'..l* r./ c )y.( g _ g o)Ip(g) 

x Jjds' f k { s' - s 0 )X P (s')e i ^ [v+Vu]i ' rk/c )e i ( fP ^ +v ‘<’’ i ')- 9P ( , ' + ‘ , >->' i 'd\ 

_ e »(#i- ^*)jfc BV /G^G7e , ( 2 ’ r l‘' +, ' , »]t i '> • r <-»o-r*l/e) 

X JJ ds e^+^-l 4 J r ‘“ p -]/«)/ (( i _ s 0 )/t(s - so)/p(s) 

= ^*'~* k) kBs/G~Gl JJ ds e <(a'[*'+«.»-*oMr.-rd/c) /<(g _ g o)/t(s _ io)/p(g) 

= &*-Mh B >/<Wh JJds e ^+^]H-io] B tk ,c) fi(i _ So)/i(g _ So)/p(g) 


(A-22) 
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i 


where the substitutions rf = s 0 r,/c, r t m = s 0 -r k /c and B,* = r* - r, have been made. 
Similarly 


x JJ ds e ‘( 2 *[*'+ ,, «»]I s - i ®l Bik / <: ) fi( s - s 0 )fk(s - so)/«/c(s) (A-23) 

In considering the effects of noise, it is assumed 10 that Eq. (A-21) can be approximated as 


N 

Pik 


^ dveW-^v? (M + T?)v* k N (v,t + t f 

k B yjT&(j* V du eK«r-*r)eff 2x ^*+^ m l +2 ' v, » T ' m+ ^ (v *‘ y '‘ +T ' m)) e _, ' (2Tw[t+Trl+2Ty, ” T ‘ m+ ^ (v ' k,, ‘ +T ‘ m)) y 
: k B\/TiTk(^J o dv e , ( 2,r, d 1 \™ -T ™l+ 2lr ‘' , “[ T *™ _T i 


- ]+ ^-_^+#w(*. lk , 1 t+Tr)-»j , («'.**.*+rr)) 




k B y/fa\\[^-^' k 


(A-24) 


where 0" is a phase randomly distributed between 0 and 2 tt, and where use has been made of the fact that 0, K,<) 
is uncorrelated between antennas and over time and frequency intervals where A uT >1. Because the phases of p jt are 
random, the double sum in Eq. (10) will be small. However, as discussed in Appendix C, noise terms become impor an 
when the process of phasing the array is considered. 


Using Eq. (A- 12) to substitute for I,/ c , 


Pii c = 2V,,lj^ d V 6{v-[uil-i^o]y (K #r) e* (#i * k) * B s/G&~k 

x JJ ds e »( 2 »I*'+‘'«»l[*- 4 °l- B ‘*]/ <! )/ j (: S - io)/fc(i - io)*(S - 

- 2 V s/c k B s/G^e i( ' H '~ Hk) e i ^ TV ‘' k ^’ /c ~ ia] Bik,c) fi(s a /c ~ s 0 )fk(s 3 / c ~ so) (A-25) 


10 This assumption is not valid if the planet or the spacecraft contributes significantly to the system temperature; in that case, the cross 
terms in Eq (A-21) will be non-negligible. However, they can be treated in a manner similar to the v v term. 
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Note that s s j c §o must be small enough that the quantity v*ly[Ss/ c ~ Sq] * / c does not change significantly over the 
interval of the time averaging. 

Similarly, 


Pik = *B\/GiG k JJ ds fi{s - s 0 )/i(s - s 0 )lp(s) jf e*( 2,r t*'+ "«o][a-S 0 ] B. k /e)^ 


(A-26) 


At this point, it is useful to define the quantity 


JJ d * ~ SO )fk(i ~ so dvS" 


i/B lfc (a-8 0 )/c) 


(A-27) 


so that 


Pit = A vkBSpy/G&leWi- 6 **)^ 


(A-28) 


?ik is a dimensionless complex quantity that depends only on the array geometry and the source structure and whose 
magnitude is always less than or equal to unity (\T? k \ -*• 1 in the short baseline limit, and \T? k \ — ► 0 in the long baseline 
limit). 

It has been implicitly assumed that the bandpass is rectangular, in which case the integral over is can be simplified: 

jf due* 2 ”** <*-*>)/•) = e >(» A'/B. i, [*-*ol/ <: )sinc ^B (i -[s - s 0 ]) (A-29) 

where sinc(ar) = sin(7ra:)/7rx. This term, often referred to as the delay beam, introduces a phase shift and lowers the 
correlation amplitude for s / so; both these effects are more pronounced for larger bandwidths A z/. 

One can therefore write 

= xk; JJ d 1 (i - i0 )/*)/,(i - s 0 )/t(s - s 0 )7p(s) sine (^B, t -[s - i 0 ]) (A-30) 


where t'tky = "la + Al//2. 

If the quantities /,, ft and sinc(At'B,* [s - s 0 ]/c) do not vary greatly over the extent of the planet, one can make 
the further simplification 

Tik = -^Urhrhuik JJ ds g) (A - 31 ) 

where f tp , /aw* are suitably weighted averages of f it f k and sinc(Ai/B, t [s - s 0 ]/c), respectively. 
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Appendix B 
Jupiter Model 


At centimeter wavelengths, Jupiter is not a simple thermal disk; there is significant synchrotron emission from the 
radiation belts and the resulting flux distribution is quite complicated [6,7]. For the purposes of this article, Jupiter 
can be modelled as the sum of three components: two circular (2-dimensional) Gaussian components (representing the 
radiation belts) and a uniform central disk. This is a simple model to work with because Eq. (A-27) can be integrated 
analytically for each of these components. 

If the brightness distribution of a source is radially symmetric about the point s c , small compared to the delay 
beam, and small compared to the primary beam of any antenna, Eq. (A-27) can be written 


?ik = 


2fifk r [sc-»o]/<0 

— /A,. k e 




(B-l) 


where r = Is - scl; /a„ is the value of the delay beam at the source; /,, /* are the average values of the antenna; 
field pattern at the source; for the *th and kth antenna; S is the total flux of the source; J(r) is the source’s radially 
symmetric brightness distribution; and Jq is a zeroth-order Bessel function. 

For a uniform disk of angular radius Rd centered at sc , 


I(r) = 


*R 2 d 


r < Rd 


= 0 r > Rp 


ik — fifkf&v ,k e 


. . cJ i (2iri/* t Bik R d / c ) 

*Kky B i* R D 


(B-2) 

(B-3) 


For a circular Gaussian of total flux S and characteristic size (1/e radius) R a , centered at s c , 


J(r) 



(B-4) 


and Eq. (B-l) becomes 


~ /i fk f A v , k e 




^^ ky B tk Ra/c) 2 


(B-5) 


Jupiter’s brightness distribution is modelled as the sum of a uniform central disk of radius Rj , centered at sj , and 
two circular Gaussians of characteristic size R B , representing the radiation belts, which are offset from the center of 
the disk by ±As B ■ The flux of the central disk is written as FdSj, where Sj is the total flux of Jupiter and Fd is the 
fraction of that flux in the disk; the integrated flux of each wing component can be written as (1 - F d )Sj/ 2. Using 
this model with Eqs. (B-3) and (B-5), 
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T- k - S- flfi. e ,( 2 ’ ri/ **v B '* [®-’-»o]/ c ) Fn c ^ 1 ( 27 rt/ »*y j 8 » t -ft J / c ) 

[ ™ C ,kyBikRj 


+ (1 — Fd) cos 




(B-6) 


If the above expression and Eq. (C-l), the condition that the array is properly phased on the spacecraft, are substituted 
into Eq. (A-28), one gets, for the cross-correlation due to Jupiter in a phased array, 

Pik = A vk B Sj VGiG k <-*»])/<) 


P cJ\ {2nu^ ky Bik Rj /c) ^ „ e ~G» c . k ,B. k R B /cf /2iri/, c *„_ 

D ’"fu.BuRj + (1 ' Fd) 7 C “ {~r L B lt As B 


In the case where s a/c = s 0 and/or v‘J ky = u e tky , this simplifies to 

Pik = AvIcbSj y/GiGt ff /j f£„ ili e < ( 2 *‘'‘ l '» B ‘ k -■•/«l/ c ) 




g i ltV iky^ikfiB / c) 2 / 27 Ti> ,C r 


+ (1 - ", - 


(B-7) 


(B-8) 


Both the total flux Sj and the flux distribution are functions of frequency and Earth-Jupiter distance. From the data 
given in [5], the model at S-band (2.3 GHz) used in this article is 


Sj 


= 6.3 


04 AU \ ' 

d ) 


Jy 


Fd = 0.3 


Rj = 24.3 


4.04 AU 
; arcsec 


|Asb| «s 2 Rj 


Rb ss 1.3f?y (B-9) 

At conjunction, d = 6.2 AU, Sj = 2.7 Jy, and Rj = 15.8 arcsec; at opposition, d = 4.2 AU, Sj = 5.8 Jy, and 
Rj = 23.4 arcsec. 

Figure B-l(a) shows a map of Jupiter from [6] at 2.9 GHz (10.4 cm), while Fig. B-l(b) shows a higher resolution 
map from [7] at 1.4 GHz. Note that in Fig. B-l(b), the disk has not been removed. Figure (B-2) shows the model of 
Jupiter used in this article in the same format as Fig. B-l(a). 
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Appendix C 
Phasing 


I. The Ideal Case 

The SNR on the spacecraft is maximized when the 
terms in the double sum in Eq. (8) add in phase. This can 
be (approximately) arranged with the appropriate choice 
of model phases </>™. By using Eq. (A-25), p*J k c will have 
zero phase for all t,Jb when 


— 4>k — 4™ — 4i + <t>T — 


27 xu\ 


•ky 


- s 0 ) 
(C-l) 


The above relation can be satisfied if, for a 11 i = l t ...N, 


The phase of the correlation function is given by 


tan$f t 


Im(m) 

Re(pjt) 


_ Im (Pjt‘) + Im(p&) + Im(pfl) 

Mp£) + R e(p5t) + Re(p%) 


(C-3) 


where Ite and Im represent, respectively, the real and 
imaginary parts of the complex function. The simplest 
phasing algorithms make the approximation that the 
spacecraft signal is the dominant contribution to i e., 


<t>? = <t>i + 


27 TV 


s/c 

sky 


®*r * (®*/c 


(C-2) 


tan *f* 




Im (Pa C ) 

R e (Pih 


(C-4) 


where r represents an arbitrarily chosen reference antenna. 

The values of <f>i are not known a priori, but there are 
N(N — l)/2 independent measurements of the phase of the 
correlation function that can be used to fit for the N model 
phases (see [8]). Note that since 6<f> iy 6<f> k include media 
effects as well as phase shifts in each antenna’s signal path, 
they vary on timescales of seconds to minutes, so that the 
model phases must be recalculated at least that frequently 
to prevent degradation of the signal. 


This approximation requires that \p\ k c \ » \pf k \ and 

\Pil C \ >> \Pil\ f° r baselines Using Eqs. (A-24), 

(A-25), and (A-28), this implies 

*V » (C-5) 

and 


II. The Effects of Random Noise and a 
Background Planet 

The above discussion ignores the contributions of ran- 
dom noise and a background planet to the observed phases. 
These complicate the process of solving for model phases, 
since what actually can be observed is the phase of the 
total correlation p ikt not the phase of p*£ c , the correlation 
due to the spacecraft signal. 


» 5\/| 


Tk j Av 
G k y T 


(C-6) 


In theory, the first of these conditions can be circumvented 
by modelling the contribution of pf k to the phase of p ik 
in the phasing algorithm. The second condition is more 
fundamental and can only be dealt with by decreasing the 
observing bandwidth (if the spacecraft is a narrowband 
source) or increasing the time over which the correlation 
is averaged. 
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Appendix D 
Weighting Factors 


The weighting factors Wi should be chosen so as to maximize /?*, which is given by Eq. (15). Such Wi will satisfy 
the condition 


= 0 (D- 1 ) 

dWj 

for all j = 1, ...N. Since one of the weighting factors can be chosen arbitrarily, this condition amounts to N — 1 
equations for N - 1 variables. 

If the effects of the planet can be ignored (i.e., S p « 0), Eq. (15) reduces to 


h 


N 


N N 


Y, w ? G *+Y,lL WiWk >/ GiGk 


i — 1 


i= 1 k?i 

Y, w ? Ti 

i= 1 


(D-2) 


and Eq. (D-l) becomes 


dds 

j^WfTi 

.1=1 

N 

2WjGj +J2w k y/GjCh 

k*j 

- 2WjTj 

’ N N N 

E wfGi + E E Wi Wt ^ G ' Gk 

i=i «=i k*i 

dWj 

E W > T ‘ 

*= i 

2 


It is not terribly difficult to show that this is satisfied for 


Wi = 


Ti 


(D-4) 


If the contribution of the planet to the system temperature at each antenna is accounted for, but the correlated 
noise terms are ignored, the optimal weighting is 


Wi = 


Vg± 

Ti + GiSp 


(D-5) 


This is appropriate for an extremely extended array. If correlated noise terms are included, it is difficult (perhaps 
impossible) to solve for Wi in closed form, but it can be done numerically. For the array configuration at the DSN 
complex in Australia, the weights obtained, including the effects of correlated noise, differ by only a few percent from 
those given by Eq. (D-5), and they have a negligible effect (<0.5 percent) on the values of ft#. In fact, though the 
values of Wi obtained, including the effects of correlated noise, often differ by ~10 percent from those obtained from 
Eq. (D-4), the resulting values of differ by <1 percent. Therefore, throughout this article, Eq. (D-4) is used to 

calculate weighting factors. 
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